Theory of polarization enhancement in epitaxial BaTiOa/SrTiOs super lattices 
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The spontaneous polarization of epitaxial BaTiOs/SrTiOs superlattices is studied as a function of 
composition using first-principles density functional theory within the local density approximation. 
With the in- plane lattice parameter fixed to that of bulk SrTiC>3, the computed superlattice polar- 
ization is enhanced above that of bulk BaTiOs for superlattices with BaTiOs fraction larger than 
40%. In contrast to their bulk paraelectric character, the SrTi03 layers are found to be tetragonal 
and polar, possessing nearly the same polarization as the BaTiOs layers. General electrostatic ar- 
guments elucidate the origin of the polarization in the SrTiC>3 layers, with important implications 
for other ferroelectric nanostructures. 



PACS: 64.70.Nd, 68.65.Cd, 77.22.Ej, 77.84.Dy 

Experimental capabilities now allow layer-by-layer epi- 
taxial growth of perovskite-based oxides, facilitating the 
exploration of a wide range of artificial materials inac- 
cessible by conventional solid-state synthesis. Recently, 
short-period BaTiOs/SrTiOs superlattices, with modu- 
lation lengths down to three perovskite layers, have been 
grown on SrTiC>3 substrates using pulsed-laser deposition 
(PLD) 1 and reactive molecular beam epitaxy (MBE). 2 
These superlattices are reportedly free of dislocations 
and coherently matched to the substrate, implying misfit 
strains of over 2% in the BaTiC-3 layers. Strains of this 
magnitude are expected to increase the polarization, 3,4 
and the superlattice geometry preserves the high-strain 
state and prevents relaxation of the BaTiC-3 (BT) layers. 
But since the BT layers are interleaved with SrTiC>3 (ST), 
which has zero spontaneous polarization in the bulk, the 
degree to which the overall polarization might be en- 
hanced in the superlattice remains in question. 

In this Letter, we address this issue directly through 
study of a series of ideal short-period BT/ST super- 
lattices with varying composition using first-principles 
density functional calculations. Our calculations show 
that these superlattices are ferroelectric, and predict 
that some possess polarizations significantly larger than 
bulk BT. An atomic-level analysis also reveals a signifi- 
cant non-zero polarization and strain in the ST layers, 
a direct result of internal polarizing fields originating 
in the BT layers. These results, consistent with recent 
measurements, 1 ' 2 illustrate the importance of electrical 
boundary conditions for sustaining polarization in these 
and other nanoscale ferroelectric materials. 

To predict the ground state structure and polarization 
in BT/ST superlattices with various numbers of BT lay- 
ers, we use density functional theory within the local den- 
sity approximation (LDA). 5 A plane- wave basis set and 
projector-augmented wave potentials 6 as implemented in 
the Vienna ab initio Simulations Package (VASP) 7,8 are 
employed. We consider period-five superlattices epitaxi- 
ally grown on ST, and construct 1x1x5 supercells, des- 
ignated 4/1, 3/2, 2/3, and 1/4, where the notation x/y 
refers to x perovskite layers of BT and y layers of ST. 



The underlying ST substrate is treated implicitly by con- 
straining the in-plane lattice constant of each superlat- 
tice to 3.863 A, the value we calculate for the equilibrium 
lattice constant of cubic ST. The computed in-plane lat- 
tice constant for bulk tetragonal BT (3.945 A) is like- 
wise slightly smaller than experiment, a well-known ar- 
tifact of the LDA. The resulting mismatch between the 
theoretical lattice constants is 2.1%, in excellent agree- 
ment with that observed experimentally. The ions within 
each supercell are allowed to relax toward equilibrium 
along [001], within space group P4mm (point group C4,,), 
until the Hellmann-Feynman forces are less than 10~ 3 
eV/A; the total energy with respect to the normal (or c- 
axis) lattice parameter of each superlattice is minimized 
concurrently. Brillouin zone (BZ) integrations are per- 
formed with a 6x6x2 Monkhorst-Pack mesh. A 44 Ry 
plane-wave cutoff is used for all calculations. For Berry- 
phase polarizations, 9 we find that 6 /c-points/string along 
[001] and 6 strings in the irreducible wedge provide well- 
converged results. 

All calculations are performed under periodic, "short- 
circuit" boundary conditions, equivalent to a metallic 
substrate with the lattice constant of ST and full charge 
compensation on top and bottom "electrode" layers. 
Canting of polarization toward [111] is found to be ener- 
getically unfavorable, although we defer to a future study 
a full investigation of the possibility of zone-boundary 
octahedral rotations, such as those present in ST at low 
temperature. 10,11 

The c/a lattice parameters computed for the four 
period-5 superlattices, as well as those for bulk ST (0/5) 
and bulk strained BT (5/0), are provided in Table I. 
With the in-plane lattice parameter fixed to the com- 
puted ST lattice constant, the c/a ratio of BT expands 
to 1.0685 (x5 layers = 5.3426 for the 5/0 superlattice), 
in agreement with the experimentally determined Pois- 
son ratio, 1 ' 2 and the c/a ratio and cell volume decrease 
by about 2% per additional SrTiOs layer. The decrease 
in c/a results in a concomitant reduction in the polar 
distortion within each layer. This is quantified by divid- 
ing each superlattice into five Ti-centered unit cells (or 
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TABLE I. Structural parameters computed for superlat- 
tices with in-plane lattice constant a = 3.863 A, and with 
1st layers of SrTi03 and Zbt layers of BaTiOs (a is their ra- 
tio). (c/a)BT,{c/a)sT, and (c/a)\ are the average local c/a 
parameters within BT, ST, and interface layers for each su- 
perlattice, as described in the text. Po=24.97 /iC/cm 2 is the 
computed value for bulk tetragonal BaTiC>3. Blank entries 
indicate layers absent for a given a. 
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layers), with A cations at the corners, and examining lo- 
cal displacements from the pseudocubic positions within 
each layer. The individual perovskite layers arc labeled 
as BT, if bounded by two BaO layers; ST, if bounded 
by SrO layers; or I, if bounded by one BaO and one SrO 
layer (one of the two interface layers within the super- 
cell). The distance between neighboring A cations can 
be regarded as a local c/a; the average values of this ra- 
tio for each layer type, (c/a), appear Table I. The large 
(c/a)eT are expected, based on the considerable in-plane 
strain. 4 More surprisingly, the average local c/a of the 
ST layers deviates slightly from unity: the ST unit cell 
is evidently expanded within the superlattice relative to 
the bulk. The strain is accompanied by relative displace- 
ments within the ST layers, where in particular we find 
the equatorial oxygens (Wyckoff positions 2c) to move 
significantly offsite. 

The observed structural trends are reflected in the 
computed superlattice polarizations, also given in Table 
I. Despite the presence of one SrO layer, the polarization 
of the 4/1 superlattice is predicted to be considerably 
enhanced over that of BT under standard conditions in 
the bulk, and also larger than expected if the ST were 
nonpolar. Replacing one BaO layer with SrO reduces the 
polarization by only about 10%, while a simple scaling 
of the polarization by the volume fraction of BT would 
result in a 20% reduction. The 3/2 and 2/3 superlat- 
tices are found to retain 80 and 64% of the polarization 
of pure strained BT, again surpassing the scaled values 
of 60% and 40%. The simple scaling argument is in- 
adequate because it assumes bulk behavior for the ST 
and strained BT layers, and neglects any electrostatic 
coupling between them. Internal polarizing fields play 
an essential role in determining the properties of these 
superlattices, as has been noted previously for other su- 
perlattice systems. 12,13 

To gain insight into the electrostatic coupling between 
layers, we decompose the total superlattice polarization 
into contributions from the individual Ti-centered per- 
ovskite layers defined above. (An alternate but equally 



FIG. 1. Local polarization enhancement (Pa/Po) by layer 
A for each superlattice. Each layer is labeled BT, ST, or I as 
discussed in the text. As before, Po=24.97 /iC/cm 2 . 

valid A cation-centered local cell choice does not result 
in significant changes in our analysis.) The local polar- 
ization P\ of cell A is expressed as a function of displace- 
ments within that cell as 

dP 1 

Px = Y.x — (u Xt -u XOt ) = -j2z;Au Xl , (i) 

where Au\i is the displacement of ion i in unit cell A, Z* 
is the dynamical effective charge of ion i, and f2 is the 
volume of the superlattice supcrcell. (All polarizations 
and displacements arc along [001].) Since there are two 
different apical O ions (corresponding to Wyckoff posi- 
tions lb) associated with the buckled AO plans bounding 
the cell, we average their contribution to each cell, i.e., 
AuA,ib=5(AwA+i,ib + AuA,ib)- A similar prescription 
has previously been used in a study of domain walls in 
PbTi03. 14 The dynamical effective charges are evaluated 
by finite differences, in a reference structure with atomic 
displacements identically zero within every layer; their 
values will appear elsewhere. 15 

In Figure 1 the local polarization is shown as a func- 
tion of layer for each superlattice. The most conspicu- 
ous feature of these profiles that the local polarization is 
nearly constant throughout each superlattice, minimizing 
electrostatic energy costs associated with the build-up of 
polarization charge V-P at the interfaces. Not only are 
the ST layers polarized, but their local polarizations are 
close to those of the BT layers. Notably, the ST po- 
larization, induced by the presence of the BT layers, is 
significant even in the 1/4 superlattice with just a single 
BaO layer, and even in the absence of in-plane strain. 
The polarization in the BT layers, on the other hand, is 
reduced and close in magnitude to that in the ST lay- 
ers. If the local polarizations of the BT and ST layers 
were taken as those of the bulk, then the electric fields 
within the constituent BT and ST layers would be finite 
and have opposite signs. The field within the BT layers 
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would oppose the polarization, and that in the ST layers 
would polarize the nominally nonpolar layers. 

To examine this in more detail, we model the system 
by two slabs of linear dielectric media, having dielectric 
constants €b and e s , in a parallel plate geometry. Short- 
circuit boundary conditions require the electric fields 
within each slab to be related by £b^BT=— £ s Ist, where 
Ibt and 1st correspond to the number of layers of BT 
and ST, as in Table I. Setting their electric displace- 
ments D = £ + AttP to be equal, and using P s = x s £ s and 
Pb = Xb£b+Pb/o> tne polarization of each of the two slabs 
can be expressed in terms of the dielectric constants of 
the two materials, the initial spontaneous polarization, 
P 5 /o, of strained BT, and their respective thicknesses. 
The average of the total slab polarization P, weighted by 
the number of layers of each constituent, then becomes 



P = 



5/0 



1 + a(e h /e s y 



(2) 



where a = 1st /Ibt and P5/0 is the computed polariza- 
tion of strained BT, here 39.20 ^C/cm 2 (as in Table I). 
This macroscopic formula, valid to linear order in the di- 
electric constants, can be seen to explain the composition 
dependence of the superlattice polarization. In Fig. 2 we 
plot the total polarization of each superlattice, calculated 
from first principles, against a. We then perform a one- 
parameter fit to this data, obtaining the dashed curve in 
Fig. 2. Notice that the curve approximates the data re- 
markably well, even though our expression for the polar- 
ization is derived from macroscopic electrostatics and the 
superlattice layers are atomically thin. Thus to maximize 
the superlattice polarization, one must increase the frac- 
tion of BT layers (while preserving the high strain state), 
consistent with the recent measurements of Shimuta et 
al. 1 It is particularly striking that even for supcrlattices 
containing as little as 40% BT, one can expect to achieve 
a polarization as large as that of bulk BT. 

The novel properties of short-period superlattice mate- 
rials are, quite generally, expected to arise from three sep- 
arate effects. The first is the sustainability of large strains 
in sufficiently thin lattice-mismatched layers. This can 
produce a significant change in properties from the bulk 
at standard conditions, manifest here in the enhanced 
polarization of the BT layers. We also expect electro- 
static effects if the two materials have different suscep- 
tibilities or a polarization mismatch (i.e., different bulk 
polarizations), as in the present case. Indeed, internal 
fields can largely determine polarization and structural 
properties, resulting in a uniform polarization through- 
out this chemically inhomogencous BT/ST system. Fi- 
nally, we anticipate changes in properties associated with 
a high concentration of interfaces, where the bonding and 
structure will in general depart from that of the interior of 
the layers or the bulk. While in this system the interfaces 
are relatively gentle, apparently playing a minor role in 
the structure and polarization, they may be important 
or even dominant in other superlattice structures, espe- 
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FIG. 2. Polarization enhancement (P/Pn), computed from 
first principles, as a function of a = l s /lb for each superlat- 
tice (filled circles). eb/e s =0.4229 provides the best fit of our 
calculations to Eq. (2). 



daily if the thickness of the individual constituent layers 
is only a few unit cells. 

In conclusion, we have shown, using a series of first- 
principles calculations, that significant polarization en- 
hancement can be achieved in perovskite oxide superlat- 
tices. This enhancement arises from the combined ef- 
fects of strain, induced in the BT layers by the epitax- 
ial growth, and internal electric fields, associated with 
the superlattice geometry, which polarize the ST layers. 
The induced polarization observed here complements a 
previous prediction of ferroelectricity in ST layers under 
epitaxial stress. 3 ' 16 Our analysis reveals important phys- 
ical factors that influence the behavior of the superlattice 
system for given constituent layers and thickness, which 
will aid in the prediction of properties of a wider class 
of systems and provide a valuable guide for the design of 
artificially-structured materials. 
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